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ABSTRACT 

We present a study of the compatibility of some current models of the diffuse Galactic continuum 7- 
rays with EGRET data. A set of regions sampling the whole sky is chosen to provide a comprehensive 
range of tests. The range of EGRET data used is extended to 100 GeV. The models are computed 
with our GALPROP cosmic-ray propagation and gamma-ray production code. We confirm that the 
"conventional model" based on the locally observed electron and nucleon spectra is inadequate, for all 
sky regions. A conventional model plus hard sources in the inner Galaxy is also inadequate, since this 
cannot explain the GeV excess away from the Galactic plane. Models with a hard electron injection 
spectrum are inconsistent with the local spectrum even considering the expected fluctuations; they 
are also inconsistent with the EGRET data above 10 GeV. 

We present a new model which fits the spectrum in all sky regions adequately. Secondary antiproton 
data were used to fix the Galactic average proton spectrum, while the electron spectrum is adjusted 
using the spectrum of diffuse emission itself. The derived electron and proton spectra are compatible 
with those measured locally considering fluctuations due to energy losses, propagation, or possibly 
details of Galactic structure. This model requires a much less dramatic variation in the electron 
spectrum than models with a hard electron injection spectrum, and moreover it fits the 7-ray spectrum 
better and to the highest EGRET energies. It gives a good representation of the latitude distribution 
of the 7-ray emission from the plane to the poles, and of the longitude distribution. We show that 
secondary positrons and electrons make an essential contribution to Galactic diffuse 7-ray emission. 
Subject headings: diffusion — cosmic rays — ISM: general — Galaxy: general — gamma rays: obser- 
vations — gamma rays: theory 
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1. INTRODUCTION 

Diffuse continuum 7-rays from the interstellar medium 
are potentially able to reveal much about the sources 
and propagation of cosmic rays (CR), but in practice the 
exploitation of this well-known connection is problem- 
atic. While the basic processes governing the CR propa- 
gation and production of diffuse 7-ray emission seem to 
be well-established, some puzzles remain. In particular, 
the spectrum of 7-rays calculated under the assumption 
that the proton and electron spectra in the Galaxy re- 
semble those measured locally reveals an excess ab ove 1 
GeV in the EGRET spectrum (J Hunter et aIlll997D . 

The Galactic diffuse continuum 7-rays are produced 
in energetic interactions of nucleons with gas via neutral 
pion production, and by electrons via inverse Compton 
(IC) scattering and bremsstrahlung. These processes are 
dominant in different parts of the spectrum, and there- 
fore if deciphered the 7-ray spectrum can provide infor- 
mation about the large-scale spectra of nucleonic and 
leptonic components of CR. In turn, having an improved 
understanding of the Galactic diffuse 7-ray emission and 
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the role of CR is essential for unveiling the spectra of 
other components of the diffuse emission and is thus of 
critical importance for the study of many topics in 7-ray 
astronomy, both Galactic and extragalactic. 

The puzzle of the "GeV excess" has lead to an at- 
tempt to re-evaluate the reaction of 7r°-production in 
pp-interactions. However, a calculation made using 
modern Monte Carlo event gen erators to si mulate high- 
energy pp-coUisions has shown ("Moril HoOTl) that the 7- 
ray flux agrees rather well with previous calculations. 
A flatter Galactic nucleon spectrum has been suggested 
as a possible solution to the "G eV excess" problem 
l|Gralewicz et al.l 119971: iMoril 11997!) . Explaining the ex- 
cess requires an ambient proton spectrum power-law in- 
dex of about -2.4-2.5, compared to - 2.75 measured lo- 
cally (for a summary of recent data see lMoskalenko et alJ 
12002(1 . Such large variations in the proton spectrum 
are, however, improbable. A sensitive test of the large- 
scale-average proton spectrum has been proposed by 
Moskalenko. Strong. & Reimer (1998) based on the fact 
that secondary antiprotons and 7-rays are produced in 
the same pp-interactions. The secondary antiprotons 
sample the proton spectrum in a large region of the 
Galaxy, and a flatter nucleon spectrum in distant regions 
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would lead to overproduction of secondary antiprotons 
and positrons. The "hard nucleon spectrum" hypothesis 
has effectively been excluded by recent measuremen ts of 
p/p ratio at high energies (e.g.. iBeach et al.ll200H) . In 
addition, new accurate measurements of the local proton 
and Helium spectrum allow less freedom for deviations 
in the Tr^-decay component. 

A "hard electron s pectrum" hypothesis has been inves- 
tigated bv IPorter fc Prothcroc (1997), Pohl & Esoosit^ 
l)1998j) . and Ahar onian fc AtovanI 1)2000) . An essen- 
tial idea of this approach is that the locally-measured 
CR spectrum of electrons is not a good constraint be- 
cause of the spatial fluctuations due to energy losses 
and the stochastic nature of the sources in space and 
time; the average interstellar electron spectrum respon- 
sible for 7-rays via IC emission (and bremsstrahlung) 
can therefore be quite different from that measured lo- 
cally. A n extensive study of this hy pothesis has been 
made bv IStrong. Moskalenko. fc Reim er (200^; in this 
model a less dramatic but essential modification of the 
proton and Helium spectrum (for the Tr'^-decay compo- 
nent) was also required. The latter was still consistent 
with the locally-observed proton spectrum, as it should 
be since the proton fluctuatio ns are expected to be small 
IjStrong fc Moskalenkall2001aj) as the result of their neg- 
ligible energy losses. The "hard electron spectrum" hy- 
pothesis suffers however from the following problems: 

• It is hardly compatible with the local electron 
spectrum even considering the fluctuations due to 
stochastic sources and ener gy losses, as shown by 
a 3D t ime-dependent studv ((Strong fc Moskalenko! 
UMa); 

• The fit to the shap e of the 7-ray spec trum is still 
poor above 1 GeV l|Strong et alil200(l : 

• It cannot reproduce the spectrum in the inner 
and outer Galaxy and intermediate /high latitudes 
simul taneously l|Strong. Moskalenko. fc Reimeri 

uni; 

These problems were already evident before the present 
study, but now we show in addition that 

• it predicts significantly higher intensities than the 
EGRET data above 10 GeV. 

Another suggestion which has been made 
IjBerezhko fc Void I2000() is that the 7-ray spectrum 
contains a contribution from accelerated particles 
confined in SNR. The SNR proton and electron spectra, 
being much harder than the interstellar CR spectra, 
produce Tr'^-decay and IC 7-rays adding to the appar- 
ently diffuse 7-rays, while the SNRs themselves are too 
distant to be resolved into individual sources. 

An alternative model involving spatial variation of 
the CR propag ation conditions has been proposed by 
lErlvkin fc mifendalel ll2002albD . 

A shortcoming of previous analyses was that the 
comparison with EGRET data was limited to particu- 
lar regions, and the rich EGRET data have remained 
not fully exploited. IHunter et al.l (^997) made an ex- 
tensive comparison of the spectra near the Galactic 
plane |6| < 10°. Other analyses have co ncentrated on 
particular molecular clouds: Ophiuchus l|Hunter et al.l 



119941) ■ Orioii llDigel et all I1999H Ceoheus and P erseus 
iPigel et alJll996H. Monoceros pigel et all 12001(1 . high 
lat itudes iS rcckumar et al."1998'). The previous analysis 
bv iStronge t al. (2000) was limited to the inner Galaxy 
at low latitudes, and profiles integrated over large re- 
gions of longitude or latitude. In that study, we com- 
pared a range of models, based on our CR propagation 
code GALPROP, with data from the Compton G amma 
Ray O bservatory. Relative to the work of Hunter et alJ 
|1997() we emphasized the connection with CR propa- 
gation theory and the importance of IC emission, and 
less to fitting to structural details of the Galactic plane. 
The study confirmed that it is rather easy to get agree- 
ment within a factor ^2 from a few MeV to 10 GeV 
with a "conventional" set of parameters, however, the 
data quality warrant considerably better fits. 

In the present paper we attempt to exploit the fact 
that the models predict quite specific behaviour for dif- 
ferent sky regions and this provides a critical test: the 
"correct" model should be consistent with the data in all 
directions. We show that a new model, with less dra- 
matic changes of electron and nucleon spectra relative 
to the "conventional" model, can well reproduce the 7- 
ray data. The changes consist in renormalization of the 
intensities of the electron and proton spectra, and a rel- 
atively small modification of the proton spectrum at low 
energies. The model is compatible with locally observed 
particle spectra considering the expected level of spatial 
fiuctuations in the Galaxy. We extend the 7-ray data 
comparisons over the entire sky and to 100 GeV in en- 
ergy. We also exploit the recent improved measurements 
of the local proton. Helium, as well as antiproton, and 
positron spectra which are used as constraints on the 
proton spectrum in distant regions. 

Our approach differs from that of ' Hunter et alJ ()1997() 
in that it is b ased on a model of CR propagation while 
IHunter et al.l use CR-gas coupling and a relativelv small 
IC component. It is also different from IStrong fc Mattoxl 
(^96) in that it is based on a physical model, while 
that work was based on model-fitting to gas surveys to 
determine the 7-ray emissivity spectrum as a function of 
Galactocentric radius. The current study concentrates 
on spectral aspects of the 7-ray emission; the question of 
the CR source gradient a nd the distribution of molecular 
hyd rogen is addressed in ((Strong. Moskalenko. fc ReimeJ 
l2004a). 

The selection of a good model for the diffuse Galac- 
tic emission is critical to another topic, the extra- 
ga lactic 7-ray backgr oun d (EGRB). We have argued 
in IStrong etaJ] ((2000(1 and iMoskalenko fc Stroll l(2000l^ 
that IC from a large halo can make up a substan- 
tial fraction of the high-latitude emission and hence re- 
duce the residual EGRB (and m odify its spectrum ve\- 
ativ e to iSreekum ar et alJ ll99Stl. In a co mpanion pa- 
per llStrong. Moskalenko. fc R,eimeill2004M we present a 
comprehensive discussion of the EGRB, with a new esti- 
mate which is used in the present paper. 

2. MODELS 

2.1. GALPROP code 

The principles of the GALPROP code for CR prop- 
agation and 7-ray emission hav e been descr i bed i n 
IStrong fc Moskalenkol l(1998() and IStrong et all l(2000(l . 



Diffuse Galactic Continuum Gamma Rays 



Table 1. Particle injection spectra and normalizations. 





ID 




Proton spect 


'um 




Electron spect 


rum 


Model 


Injection 
index^ 


Break 
rigidity, GV 


Normalization'' 
@ 100 GeV 


Injection 
index'' 


Break 
rigidity, GV 


Normalization" 
@ 32.6 GeV 


Conventional 
Hard electron 
Optimized 


44.500180 
44.500181 
44.500190 


1.98/2.42 
1.98/2.42 
1.50/2.42 


9 
9 
10 


5.0 X 10-2 
5.0 X 10-2 
9.0 X 10-2 


1.60/2.54 

1.90 
1.50/2.42 


4 
20 


4.86 X 10-3 
1.23 X 10-2 
2.39 X 10-2 



Note. — The GALPROP model IDs are given for future reference; the corresponding parameter files contain a complete specification 
of the models. 
^Below/above the break rigidity. 
'^ Normalization of the local spectrum (propagated). Units are m-2 sr-^ s~^ GeV-^. 



Since then the co de^ has bee n entirely re-written in 
C++ ( Moskalcnko et aP 120021 and references therein) 
using the experience gained from the original (fortran) 
version, with improvements in particular in the genera- 
tion of 7-ray skymaps. Both 2D (radially symmetric) 
and full 3D options are available, the latter allowing 
also ex plicit time-dependence with st ochastic SNR source 
events f St rong fc Moskalenkoll2001bD . For this paper the 
2D option is sufficient since we need only kpc-scale av- 
eraged CR spectra (even if these differ from local CR 
measurements) . 

An important point to note is that even in the 2D case, 
the symmetry applies only to the CR distribution; for the 
gas-related components (7r°-decay and bremsstrahlung) 
of the 7-ray skymaps we use 21-cm line survey data for 
H I and CO ( J = 1 ^ 0) survey data for H2, in the form 
of column densities for Galactocentric "rings," using ve- 
locity information and a rotation curve (see Appendix 
for details) . In this way details of Galactic structure are 
included at least for the gas, at a sufficient level for the 
present limited state of knowledge on e.g. the relation 
of cosmic rays to spiral structure. The longitude range 
350° < Z < 10° is not included in the H i and CO sur- 
vey data due to lack of kinematic information; for the 
analysis interpolated values are used, and this is found 
to be fully consistent will the 7-ray data. The interstel- 
lar radiation field (ISRF) for computing IC emission and 
electron energy losses is the same as that described and 
used in Strong ct al. (2000); pending a new calculation 
(an ambitious project) this is the best we have available. 
Although the uncertainty in the ISRF is a shortcoming, 
note that since we fit to the 7-ray data by adjusting the 
electron spectrum, inaccuracies in the ISRF spectrum 
will tend to be compensated. 

Th e radial distribution of CR sources used is the same 
as in lStrong et al.l l|2000() . since we find this empirically- 
derived form still gives a good reproduction of the 7- 
ray longitude distribution"^. Although flatter than the 
SNR distribution (e.g., Case & Bhattacharva 1998), this 
may be compensated by the gradient in the CO-to-Il2 
conversion factor whose metallicity and temperature de- 
pendences have the net effect of causing the factor to 
increase with R IjPapadopoulos. Thi. fc Vit]|l2002t llsraell 

2 As usual the code and documentation is available at 

Etpi77www.mpe.mpg.de/~aws/aws.html 

•^ I'b r earlier work on the CR distribution seeJStecker^Jone^ 
ilT977l). liTar ding fc Stecker (1985), Bloe men et al.l J198fiiK 
■Strong et alT>198ai) . lStrong fc Mattoxi >199fli') . 



119971) . We use a uniform value of Xr.n= 1- 9 x10^ 
mol ecules cm-^/(K km s~^ ) as in iStrong et al.l l|200(]n 
and "S trong fc Mattoja ((1996^ ; this is consistent with the 
value (1.8 ± 0.3) X 10^° molecules cm-^/(K k m s^^) from 
a rece nt (non-7-ray) CO survey analysis bv lDame et alJ 
ll200lh . 




Kinetic energy, GeV/nucleon 

Fig. 1. — B/C ratio as calculated in reacceleration 

model. Lower curve - LIS, upper - modula ted (j> = 

450 MV). Data below 200 MeV/nucleon: AC E IDavis et all 
l200(1) . Ulysses CDuVernois. Simpson, fc Thaverf 119961) . Voy- 
ager (Uukasiak. McDonald. & Webber 199(^i high energy data: 
HEAO-3 ^E ngclmannctaLi tlflSI^), for other references see 
■Stephens fc Streitmatterf 119981) . 

The parameters of the models are summarized in Ta- 
ble n The models differ only in the injection spec- 
tra of protons (and He) and electrons, while the injec- 
tion spectra of heavier nuclei are assumed to have the 
same power-law in rigidity, for all models. For propa- 
gation, we use essentially the same di ffusion reaccelera- 
tion model, model DR, as described in'Mosk alenko et alJ 
(|2002) . The propagation parameters have been tuned to 
fit the B/C ratio (F ig. [TJ using improved cross-sections 
IjMoskalenko fc Mashnikii2003^1 . The spatial diffusion co- 
efficient is taken as PDo{p/po)^, where Dq = 5.8 x 10^* 
cm s""'^ at po = 4 GV, and S = 1/3 (Kolmogorov spec- 
trum). The Alfven speed is va = 30 k m s^-^. The 
halo height is taken as z^ = 4 kpc as in IStrong et alJ 
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Table 2. 


Sky regions used 


FOR COMPARISON OF MODELS WITH DATA. 


Table 3. 


Comparison 


OF MODELS USING 


x' 


FOR 


FULL SKY. 


Region 
H 


I, deg 

300-60 

330-30 

30-330 

90-270 

0-360 

0-360 

0-360 


|6|,dcg 

0-10 

0-5 

0-5 

0-10 

10-20 

20-60 

60-90 


Description 
Hunter et al. (1997) reei 


3n 

nner Galaxy 


Energy range 
MeV 


44.500180 
conventional 


44.500181 
hard electron 


44.500190 
optimized 


Number of 
data points 


A 
B 
C 
D 

E 
F 


inner Galaxy 
Galactic plane avoiding i 
outer Galaxy 
intermediate latitudes 1 
intermediate latitudes 2 
Galactic poles 


30-50 

50-70 

70-100 

100-150 

150-300 

300-500 

500-1000 

1000-2000 

2000-4000 

4000-10000 

10000-20000 

20000-50000 

30-50000 


90 

19 

18 

38 

33 

43 
140 
382 
441 
247 

56 

22 

1528 


38 
31 
49 
89 
82 
24 
59 
216 
243 
53 
54 
35 

974 




34 
26 
30 
47 
41 
24 
32 
61 
93 
49 
21 
4 

462 




78 
78 
78 
78 
78 
78 
78 
78 




in accordance 
/primary ratios ( 


with our analysis of 


CR sec- 


76 
66 
23 

7 


(laooo), 




ondary 


Moskalenko. Mashnik. 


& Strond 


796 


120011 Moskalenko et al 


120021 and references 


therein) . 





However, values Zh differing by 50% (the estimated er- 
ror) with corresponding adjustment of Dq would provide 
essentially similar results since the IC contribution scales 
mainly with the electron spectrum which is here treated 
as a free parameter. 

The spectra are compared in the regions summarized 
in Table 121 Region A corresponds to the "inner radian", 
region B is the Galactic plane excluding the inner radian, 
region C is the "outer Galaxy", regions D and E cover 
higher latitudes at all longitudes, r egion F is "Galacti c 
poles" . Region H is the same as in iHunter et al." (1997) 
and is used for comparison with results of Hunter et al. In 
addition to spectra, profiles in longitude and latitude are 
an essential diagnostic; our latitude profiles are plotted 
logarithmically because of the large dynamic range from 
the Galactic plane to the poles. 

The EGRB used her e is bas ed on the new determi- 
nation by Str ong et al.l lj2004bj) . Since this was derived 
for the EGRET energy bands, it is interpolated in order 
to produce a continuous spectrum for combining with 
the model Galactic components. The present analysis is 
however not sensitive to the details of the EGRB. Since 
our COMPTEL data do not contain the EGRB (see Sec- 
tion I^J, we do not extrapolate the EGRB beyond the 
EGRET energy range when comparing with data. 

2.2. Presentation of results 

The output of the GALPROP runs is in the form of 
FITS files; the visualization^ in the form of spectra and 
profiles, and comparison of the results with data involves 
integrations over sky regions and energy as well as con- 
volution. The predicted model skymaps are convolved 
with the EGRET p oint-spread function as described in 
IStrong et al.l IJ2000D . For the profiles the convolved model 
is directly compared with the observed intensities. For 
the spectra, the procedure is slightly different: the pre- 
dicted (unconvolved) intensities are compared with in- 
tensities corrected for the effect of convolution as given 
by the model under study. This procedure has the advan- 
tage that the spectra are spatially deconvolved, allowing 
for more direct interpretation and also the combination of 
data with other experiments, such as COMPTEL, with 
different instrument response functions. The effect of 

"^ An additional program GALPLOT has been developed for this 
purpose, with flexible plotting options and convolution; this will be 
made available with future versions of GALPROP. 



this procedure on the spectra is only significant below 
500 MeV. 

2.3. Statistical test 

The choice of model in this work is mainly subjective, 
based on visual inspection of spectra and profiles. In 
order to give also an objective criterion, a x^ statistic 
has been computed over the full sky, for each EGRET 
energy range between 30 MeV and 50 GeV (TablcEl ); The 
binn ing for this test is the same as used in .Strong et alJ 
|2Q.Q43) : raster scanned bins in longitude for latitud e 
width 2°, giving 78 sky bins. As in IStrong etlll l)2004bj) . 
a lower limit of 10 counts per bin were accepted, so at 
high energies the number of bins is reduced. The error 
is computed as the sum of the statistical error and the 
systematic error as described in section IXTI 

3. 7-ray AND COSMIC RAY MEASUREMENTS 

3.1. EGRET data 

We use the co-added and point-source removed 
EGRET counts and exposure maps in Galactic coordi- 
nates with 0.5° X 0.5° binsize at energies between 3 
MeV and 10 GeV, as described in IStrong et "all l)2000|) . 
Apart from the most intense sources, the removal of 
sources has little influence on the comparison with mod- 
els. For the spectra, the statistical errors on the EGRET 
data points are very small since the regions chosen 
have large solid angle; the systematic error dominates 
and we have conservatively ado pted a range J=15% in 
plotting the o bserved spectra (|Sreekumar et al.1 119981: 
lEsDosito et al.l ll999). For longitude and latitude pro- 
files only the statistical errors are plotted. In addition 
we use EGRET data in the energy ranges 10-20, 20-50 
and 50-120 GeV. Because the instrumental response of 
EGRET determined at energies above 10 GeV is less cer- 
tain compared to energies below 10 GeV, it is necessary 
to account for additional uncertainties. In particular the 
EGRET effective area can only be deduced by extrapo- 
lation from the calibrate d effective area at lower energies 
({Thompson et al.lll993aj) . We accordingly adopt values 
of 0.9, 0.8, and 0.7 times the 4-10 GeV effective area, 
respectively. On top of the statistical and systematic 
uncertainties as described above we account for the un- 
certainties due to the uncalibrated effective area of the 
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EGRET telescope above 10 GeV with an additional sys- 
tematic error of ±5%. However, the actual number of 
photons >10 GeV is small: 1091, 362 and 53 events re- 
spectively, and concentrated mainly in the inner Galaxy; 
hence the comparison with models above 10 GeV can 
only be made in this region. 

At low energies the EGR ET effective area include s the 
so-called "Kniffen factor" l|ThomDSon et al.lll993b(l de- 
rived by fitting the Crab spectrum; this additional un- 
certainty (factor = 2-3.4 for 30-50 MeV and 1.2-1.6 for 
50-70 MeV) should be borne in mind when comparing 
models with EGRET data. 

3.2. COMPTEL data 

The intensities are based on COMPTEL maximum en- 
tropy all-sky maps in the en ergy ranges 1-3, 3- 10 and 
10-30 MeV, as published in iStrong et"an l(T999|) . The 
intensities are averaged over the appropriate sky regions, 
with high latitudes being used to define the zero level. 
COMPTEL data is only used for the inner Galaxy spec- 
tra, since the skymaps do not show significant diffuse 
emission elsewhere. For this reason, the COMPTEL data 
shown in the figures does not include the EGRB. 

3.3. Cosmic rays 

In our conventional model we use the locally-observed 
proton. He, and electron spectra (solid lines in Figs.|21ISl). 
The nucleon data are now more precise than those which 
were available for IStrong et alJ |2000.) . The proton (and 
Helium) injection spectra and the propagation parame- 
ters are chosen to reproduce the most recent measure- 
ments o f primary and secondary nuclei, as described in 
detail in lMoskalenko et al. (2002) . The error on the dom- 
inant proton spectrum in the critical (for T r'^-decay) 10- 
100 GeV range is now only -- 5% for BESS IjSanuki et all 



|2003). Relative to protons, the contribution of He in CR 
to the 7-ray flux is about 17%, and the CNO nuclei 
in CR contribute about 3%. The He nuclei in the ISM 
contribute about 25% relative to hydrogen for the given 
ratio He/H = 0.11 by number. The total contribution of 
nuclei with Z > 1 is about 50% relative to protons. 

In our optimized model we use the proton and He high- 
energy spectral shape derived from the local data (dotted 
lines in Figs. 13 El. We allow however for some devia- 
tions in t he normalizat ion. The antiproton IjOrito et alJ 
2000; Bea ch et al.l200lD and positron data provide an im- 
portant constraint fMoskale nko et alJIlQQa: IStrong et alJ 
2P00) on the proton spectrum on a large scale. Since 
the low-energy protons and nuclei are undetectable in 
the ISM, we allow more freedom in the proton and He 
spectrum below 10 GeV. We introduce a break at 10 
GeV which enables us to fit the 7-ray spectrum while 
still remaining within the constraints provided by the 
locally-observed antiproton and positron spectra. The 
deviations from the local measurements at low energies 
can be caused by the effect of energy losses and spa- 
tial fluctuations in the Galaxy. The modification of the 
low-energy proton spectrum may also be partly a com- 
pensation for errors in the m odels of neu tral pion pro- 
duction at low energies (e .g.. lSte ckeij|19 70^) . which rely 
on the data of 1960's fsee lDermedll986l and references 
therein) and do not provide the required accuracy now. 
Besides, the low-energy protons are strongly affected by 
solar modulation; while the effect of solar modulation 
is not fully understood, it is essential below 10 GeV. 
We refer a reader to Section |H1 where various aspects 
of the uncertainties are discussed in more detail. For 
electrons, the injection index nea r ~1.8 at ~1 GeV is 
consistent (see Str ong et al.1 '2000) with observations of 
the synchrotron in dex (3 = 2.40 - 2.55 for 22-408 MHz 
lIRoger et al.lll999D and (3 = 2.57 ± 0.03 for 10-100 MHz 
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IjWebber. Simpson, fc Canelll980(l . 

Secondary an d tertiary antipro t ons a re calculated 
as described in iMoskalenko et al.1 l)2002|) . Secondary 
positron and electron pr oduction is computed usin g 
the formalism described in "Mos kalenko fc Strond l)1998h . 
which includes a reevaluation of the secondary tt - and 
X^'^-meson decay calculations. Antiprotons, positrons, 
and electrons including secondary electrons are propa- 
gated in the same model as other CR species. 

4. CONVENTIONAL MODEL 

We start by repeating the test of the "conventional" 
model; the 7-ray spectra in the 7 test regions are shown 
in Fig. 01 As required by the "conventional" tag, the 
proton and electron spectra are consistent with the lo- 
cally observed spectra (F igs. El Et- This is the same 
"conventional" model as in lStrong et al. (2000), with up- 
dated nucleon spectra, but because we co mpare with a 
more complete set of EGRET data than in lStrong et al.l 
(|2PG0) , the discrepancies become more explicit, and we 
can check whether they arise only in particular sky re- 
gions. Note that IC plays only a minor role in this type 
of model. As found in previous work, the GeV energy 
range shows an excess relative to that predicted; what 
is now evident is that this excess appears in all lati- 
tudes /l ongitude ran ges. This is consistent w i th the re- 
sults of lHunter et~a l. (1997| and iDJgel etlll l|200H) . It 
already shows that the GeV excess is not a feature re- 
stricted to the Galactic ridge or the gas-related emission. 
Further it is clear that a simple upward rescaling of the 
TT^-decay component will not improve the fit in any re- 
gion, since the observed peak is at higher energies than 
the 7r°-decay peak. In other words, since the spectrum 
is very different from Tr^-decay even at intermediate lat- 
itudes, a substantial IC component is required. The x^ 
values (Table O confirm the visual conclusion that this 
model is unacceptable. 

Note that this version of the "conventional" spec- 
trum is nevertheless in rather be t ter a greement with 
EGRET data than in ' Strong etal\ l)2000D . due to inclu- 
sion of secondary positrons/electrons, general improve- 
ments in the model (e.g., Tr^-decay, improved gas data) 
and the EGRET data treatment (direct use of the count 
an d exposure data instead of the model-fitting analysis 
of iStrong fc Mattoxl 1199 6^. The improvement is espe- 
cially evident in the 30-100 McV range, where secondary 
positrons/electrons make a substantial contribution (see 
Section O. 

A test against antiproton and positron data also shows 
"excesses" . The conventional model with reacceleration 
is known (Moskalenko et al. 2002) to produce a factor of 
^^1.5 (~2. 5g) less a ntiprotons at 2 GeV than measured 
by BESS IjOrito et aL .20003 . The antiproton spectrum 
for the conventional model is shown in Fig. O Positron 
data, though scattered, also show some "excess" at high 
energies (Fig. jS}. It is thus clear that the "excesses" in 
GeV 7-rays in all directions, in GeV antiprotons, and in 
positrons above several GeV found in the conventional 
model indicate that the average high-energy proton flux 
in the Galaxy should be more intense or our reacceler- 
ation model is invalid or there is a contribution from 
unconventional sources (e.g. dark matter). For more 
discussion of antiproton and positron t ests see Section (71 

In the "SNR source" scenario of iBerezhko fc V61M 



(|2000(1 the 7-ray spectrum in the inner Galaxy is at- 
tributed to an additional population of unresolved SNR, 
but this component cannot explain the excess observed 
at high latitudes, and hardly in the outer Galaxy^. The 
presence of the GeV excess in a ll sky regions is also 
aproblem for the suggestion by lAharonian fc Atovanl 
(j^OO) of a hard proton spectrum in the inner Galaxy. 
These explanations are therefore by themselves insuffi- 
cient, although they could give a contribution. 

5. HARD ELECTRON INJECTION SPECTRUM MODEL 



This model is essentially as in lStrong et alJ l)2000j) . re- 
computed with the current GALPROP code. The main 
feature is the electron injection index of 1.9. Compari- 
son of the spectra in the 7 sky areas (Fig. O show that 
this model reproduces the GeV excess except in the inner 
Galaxy (region A) where it is still too low. However the 
spectral shape is not well reproduced. More significant, 
comparing with the new EGRET data above 10 GeV in 
the inner Galaxy, the spectrum is much too hard. The 
X^ values (Table OJ confirm the visual conclusion that 
this model is only marginally unacceptable. 

Fig.jJlcompares the locally observed electron spectrum 
with that from the model; the deviation at high ener- 
gies is much larger than expe cted from the 3D study by 
IStrong fc Moskalenkol l)2001bfl . As discussed in the In- 
troduction, there are therefore a number of reasons to 
lead us to consider this model as after all untenable. 

6. OPTIMIZED MODEL 

Since the conventional model fails to reproduce the ob- 
served 7-ray spectrum, and the hard electron spectrum 
model is untenable, we use the diffuse 7-rays themselves 
to obtain an optimized solution. The average interstel- 
lar electron spectrum is sufficiently uncertain that we 
can look for a "solution" involving a less drastic change 
in the electron injection spectrum than the hard elec- 
tron injection spectrum model. We find that an injection 
spectrum of electrons with a steepening from 1.5 to 2.42 
at 20 GeV (see Table P) produces sufficient curvature in 
the inverse Compton spectrum to explain the observed 
shape of the 7-ray spectrum, provided the electron spec- 
trum is suitable normalized upwards by a factor of about 
4 relative to the locally observed spectrum,. The proton 
injection spectrum is also normalized upward, by a fac- 
tor 1.8; it has the same shape as for the electrons, as 
a function of rigidity, but the break energy is 10 GeV 
instead of 20 GeV. It has exactly the same slope above 
10 GeV as the conventional proton spectrum. (The pro- 
ton re- normalization factor 1.8 is not taken ad hoc but 
is chosen to reproduce the antiproton data, see SectionO 
for more details). 

The 7-ray spectra in the 7 test regions are shown in 
Fig. (HI The fits to the observed 7-ray spectra are bet- 
ter than for the conventional and hard electron spectrum 
models, both in the 1-10 GeV region and above 20 GeV. 
The spectra in different regions are satisfactorily repro- 
duced and there is no longer a significant GeV excess. 
Hence the spectrum can now be reproduced from 30 MeV 
to 100 GeV. The proposed scenario implies a substantial 
contribution from IC at all energies, but especially be- 



^ IBerezhko fc Volkl 1200(1) did not address the question of regions 
away from the inner Galaxy. 
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low too MeV and above 1 GeV. Also fC dominates at 
latitudes \b\ > 10° at all energies. 

Longitude profiles at low latitudes are shown in Fig.^ 
The agreement with the EGRET data is generally good 
considering that the model does not attempt to include 
details of Galactic structure (e.g., spiral arms), and the 
systematic deviations reflect the lack of an exact fit to 
the spectra in Fig.|Sl The largest deviation (^20%) is at 
2-4 GeV, but this is still compatible with the systematic 
errors of the EGRET data. Latitude profiles in the lon- 
gitude ranges 330° < ^ < 30°, 30° < / < 330° are shown 
in Figs. ^1 E] where the logarithmic scale is chosen 
given the large dynamic range and to facilitate the com- 



parison at high Galactic latitudes. The agreement with 
EGRET is again good, in particular the reproduction of 
the high-latitude variation confirms the importance of 
the IC component which is much broader than the gas- 
related Tr'^-decay and bremsstrahlung emission. In the 
inner Galaxy (Fig. II 0|) there is evidence for an excess 
at intermediate latitudes, perhaps related to an underes- 
timate of the interstellar radiation field in the Galactic 
halo, or special conditions in the Gould's Belt. The outer 
Galaxy latitude profiles (Fig. Ill|l are in excellent agree- 
ment with the data. 

The x^ values (Table |3J confirm the visual conclusion 
of the improvement of this model over the conventional 
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and hard-electron spectrum models. 

The local electron spectrum (Fig.O is compatible with 
the direct measurement considering fluctuations due to 
energy losses and stochastic sources and propagation 
IjStrong fc Moskalen ko 2001b), and in addition uncer- 
tainties in solar modulation at low energies. In fact the 
agreement can be even better if we consider the uncer- 
tainty in the ISRF, which can well be a factor 2 higher 
than our estimate. The electron spectrum is consistent 
with the synchrotron spectral index data IjStrong et al.l 
120001) ■ since it differs from the conventional model essen- 
tially only in the normalization, and this is in turn con- 
sistent with synchrotron. The interstellar proton spec- 
trum (Fig. ^ is also compatible with direct measure- 
ments; the factor 1.8 may be attributed to fluctuations 
over the Galaxy relative to the local value, and also to 
the uncertainty in the large-scale CR gradient. 

Below 30 MeV the predicted spectrum lies about a 
factor 2 below the C OMPTEL data, as found previously 
l|Strong et al."200Q') . There we proposed that a contribu- 
tion from compact sources is the most lik ely explanation . 
Recent results from INTEGRAL ( Stro ng et al.ll2008b(l 
indeed indicate a large contribution from sources in the 
hard X-ray band, and this would be consistent with the 
MeV region marking a transition from source-dominated 
to diffuse-dominated ridge emission. 

7. SECONDARY ANTIPROTONS, POSITRONS, 
ELECTRONS 

7.1. Tests oj the nucleon spectrum 

A sensitive test of the proton spectrum using t he p/p 
ratio has been proposed IjMoskalenko et all 1 19981) based 
on the fact that secondary antiprotons and 7-rays are 
produced in the same pp-interactions. Positrons are 
also produced in pp-coUisions and thus may be used to 
support the conclusions made from the antiproton test. 
While some deviation from the locally-observed spectrum 



of primary protons is acceptable, secondary antiprotons 
(and partly positrons) trace the primary proton spec- 
trum on scales up to ~10 kpc over the Galaxy, and hence 
allow us to put limits on deviations from the local mea- 
surements. 

Antiprotons and positrons were originally used to ex- 
clude the possibility of a har d proton spectrurn as the 
origin of the 7-ray GeV excess llStrong et al.l2000|) . How- 
ever, even for a conventional nucleon spectrum, a prob- 
lem appears in the reacceleration model in the simul- 
taneous fitting of secondary/primary nuclei ratios and 
the antiproton spectrum; the former fixes the propaga- 
tion parameters which can be used to predict antipro- 
tons, but using the measured proton spectrum leads 
to too few antiprotons. To reproduce both the sec- 
ondar y/primary nuclei r a tios a nd the antiproton spec- 
trum, iMoskalenko et al.l l)2002D sugg ested a change in 
the propagation mode at low energies. IMoskalenko et all 
f2003) proposed a contribution to the primary CR spec- 
trum from the "local bubble" as a possible solution. Our 
present model provides an alternative to th ese so lutions. 

As was not ed in iStrong et alJ l)2000|) and 
IMoskalenko et all l)2003l and references therein), 
the "GeV excess" in 7-rays and underproduction of 
antiprotons in the reacceleration model may indicate 
that the nucleon spectrum typical of large regions of 
the Galaxy differs moderately from the local measure- 
ments. The problem with secondary antiprotons in 
th e reacceleration model has been extensively discussed 
in IMoskalenko et all l)2002l I2003J) . It is apparent that 
if the solution of the 7-ray GeV excess can not be 
found in modifications of the electron spectrum alone, 
the required modifications in the nucleon spectrum 
must satisfy the constraints from both antiprotons and 
positrons. 

Figs. [SI and ini show the antiproton and positron fluxes 
as calculated in the conventional and optimized models. 
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Tfie modifications of the nucleon spectrum introduced 
in the optimized model appear to be exactly what is re- 
quired to reproduce both, antiproton and diffuse 7-ray 
data, and the positron spectrum also agrees at high en- 
ergies. At low energies the calculated positron spectrum 
is rather high but the solar modulation is a factor of 
^1000 at these energies, and besides the scatter in the 
positron data may indicate large systematic errors. 

7.2. Gamma-rays from, secondary positrons and 
electrons 

Secondary positrons in CR produced in interactions of 
energetic nucleons with interstellar gas are usually con- 



sidered a minor component of CR. This is indeed so in the 
heliosphere where the positron to all-lepton ratio is small 
at all energies, e'^ /ctot ^ 0.1. However, the secondary 
positron flux in the interstellar medium is comparable 
to the electron flux at relatively low energies '--^1 GeV 
because of the steeper spectrum of secondary positrons. 
The spectrum of secondary positrons and electrons de- 
pends only on the ambient spectrum of nucleons and the 
adopted propagation model. Figs. 13 El show the spectra 
of electrons and secondary positrons for the conventional 
and optimized models. Secondary positrons contribute 
about half of the total lepton flux at ~1 GeV. Secondary 
electrons add up another 10% fFig. I12|l . This leads to 
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a considerable contribution of secondary positrons and 
electrons to the diffuse 7-ray flux via IC scattering and 
bremsstrahlung and significantly increases the flux of dif- 
fuse Galactic 7-rays in MeV range. Therefore, secondary 
positrons (and electrons) in CR can be directly traced in 
7-rays! 

Fig. El shows the contribution of secondary positrons 
and electrons to the IC emission and bremsstrahlung. 
Secondaries contribute more than 20% of the total IC in 
the 1-10 MeV energy range. More dramatic is the case 
of bremsstrahlung, where secondaries contribute about 
60% of the total below ~200 MeV. It is the contribution 
of secondaries which improves the agreement with data 



below some 100 MeV. 

However the secondaries are not sufficient to explain 
the excess in the 1-30 MeV range observed by COMP- 
TEL, so that an additional point-so urce contribution t o 
the emission is still required here IjStrong et al.l 12000"). 
Evidence for such a point source co ntribution has indee d 
recently been found by INTEGRAL IjLebrun et al.l2004(l . 

8. DISCUSSION OF UNCERTAINTIES 

We do not discuss here possible calculation errors. 
Derivation of such errors is a very complicated matter 
given the many uncertainties in the input. Those most 
probable are the uncertainties in the Tr'^-production in 
pp-collisions at relatively low energies, nuclear cross sec- 
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Fig. 12. — Electron spectra for conventional (solid) and opti- 
mized models (dots). Upper curve - LIS, lower - modulated to 
600 MV. Secondary electrons arc shown separately for the same 
models. 
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Fig. 13. — 7-ray spectrum of optimized model with (thick lines) 
and without (thin lines) primary electrons, to show the contribu- 
tion of secondary electrons and positrons. Brtot and Br2 labels de- 
note the total bremsstrahlung and the separate contribution from 
secondary leptons, correspondingly. Similarly, ICtot, IC2 indicate 
the total IC and the contribution from secondaries. 



tions, gas distribution in the Galaxy, ISRF, systematic 
errors in the CR measurements, hehospheric modulation, 
etc. Some possible errors and their effects have been dis- 
cussed in (Moskalenko et al. (2002, 2003). Here we qual- 
itatively mention what we think may affect our conclu- 
si ons and wha t may not. 

iMoril l)1997l and references therein) has re-evaluated 
the Tr'^-production in pp-collisions using modern Monte 
Carlo event generators HADRIN, PYTHIA, and 
FRITIOF. The HADRIN code, which is designed to re- 
produce nuclear collisions at laboratory energies below 
5 GeV and describes the threshold and resonance be- 
havior of inelastic hadron-nucleon interactions, shows 
good agreement with isobar model calculations (Stecker 
Il970|) at proton kinetic energy Tp = 0.97 GeV. The iso- 
bar model is shown to reproduce the data on the sec- 
ondary 7r°-production at low e nergies, in part icular, at 
T^ = 0.56,0.65,0.97,2.0 GeV lIDermeri Il986l and ref- 
erences therein). However, the data on Tr'^-production 
at GeV energies are now 40 years old, they have large 
statistical errors and are very scattered indicating pos- 
sibly large systematic uncertainties. Given the lack of 
new data, the deviations from the isobar model calcu- 
lations by a factor of ~2 would be also consistent with 
the old data. However, the comparison to the Galactic 
diffuse 7-ray emission is now rather precise, and the un- 
certainty in the Tr'^-production at low energies may be 
critical. Our required flattening of the proton spectrum 
below 10 GeV could thus be understood as a compen- 
sation for errors in the Tr'^-prod uction physics. At high 
energi es, a comparison of Badhwar. Ste phens, fc GoldenI 
l)1977() and IStephens & Badhwar (1981) scahng model 
with resuls of PYTHIA and FRITIOF shows generally 
a good agreement, but all of them overpredict the cross 
sections at high rapidities. Though it may result in sys- 
tematic uncertainties of the 7-ray flux above ~100 GeV, 
this, however, is of less concern given the large error bars 



in the EGRET data in this energy range. 

Possible errors in the cross section of the CR nuclei 
affect the derived propagation parameters such the diffu- 
sion coefficient, Alfven speed etc. While they may be im- 
portant for calculation of CR isotopic abundances, they 
do not affect much the calculation of the diffuse 7-rays 
as we normalize the particle spectra to the given local 
values. 

Errors in the Galactic gas distribution are not so im- 
portant in the case of stable and long-lived nuclei. Such 
errors are compensated simultaneously for all species by 
the corresponding adjustment of the propagation param- 
eters (diffusion coefficient). In case of 7-rays, we compare 
with the large sky regions so that the error in the col- 
umn density in any particular direction produces a minor 
effect. 

For the calculation of the spectrum of 7-rays arising 
from IC scattering and electron energy losses, the full 
ISRF as function of position and wavelength is required, 
which was n ot available in the l iterature. The ISRF was 
evaluated in lStrong et alJ l)2000() using emissivities based 
on stellar populations and dust emission. Given the com- 
plicated and uncertain input in this calculation a factor of 
two error is quite possible. The inaccuracies in the ISRF 
are compensated in our model by adjustment of the CR 
electron flux. Therefore, if the ISRF energy density is in 
reality higher, it will result in lower normalization of the 
electron injection spectrum making it closer to the local 
one. 

Our knowledge of the heliospheric modulation is still 
incomplete, and it remains the source of a large uncer- 
tainty in the propagation models. Over the last years 
Ulysses made its measurements at different heliolatitudes 
so we know more about the solar magnetic field configu- 
ration and the solar wind velocity distribution. However 
the modulation parameters are usually still determined 
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based on the assumed ad hoc interstellar nucleon spec- 
trum. Estimates of the modulation made using the sim- 
plest force-field approximation show that the modulation 
changes the proton intensity below 1 GeV by a factor of 
10 (Fig. EJ, and by a factor of 100-1000 in case of elec- 
trons and positrons (Figs. 13 EJ- This makes it difficult 
to speculate about the reasons for deviations of the cal- 
culated spectrum from the measured one by a factor of a 
few at low energies. Modulation is important at energies 
below ~10 GeV/nucleon, while it is negligible at higher 
energies. The 7-rays with energy >10 GeV are produced 
by protons of >100 GeV and by electrons of >10 GeV 
where modulation has (almost) no effect. Lower energy 
7-rays may be affected by the uncertainties in the solar 
modulation, but this is compensated by the adjustment 
of the interstellar spectra. In our optimized model, the 
diffuse 7-rays themselves are used to constrain the in- 
terstellar particle spectra at low energies, while the con- 
straints from the local measurements are relaxed. 

Finally, due to the random nature of SN explosions 
the CR sp ectrum fluctua tes in space and time (see simu- 
lations in iStrong fc Mo skalcnko 2003)- Further, since 
more CR sources are concentrat ed in the spiral arms 
(e.g., iCase fc Bhattacharvalll99q) . the CR intensity in 
the arms might be higher, while the sun is located in the 
interarm region. These effects can cause the locally mea- 
sured CR intensity to differ from the large-scale average 
by the required factor of 2-4. 

9. CONCLUSIONS 

We have revisited the compatibility of diffuse Galac- 
tic continuum 7-ray models with the EGRET data. We 
confirm that the "conventional model" based on the lo- 
cally observed electron spectrum is inadequate, for all 
sky regions. A conventional model plus hard sources in 
the inner Galaxy is also inadequate, since this cannot 
explain the excess outside of the Galactic plane. Models 
with a hard electron injection spectrum, while reproduc- 
ing the EGRET spectrum in the few GeV region over 
much of the sky, are not compatible with the locally ob- 
served electron spectrum (the expected fluctuations are 



not sufficient) and are inconsistent with EGRET data 
above 10 GeV. 

A new model, with relatively mild deviations of the 
electron and proton spectra from local, is shown give a 
good reproduction of the diffuse 7-ray sky. The agree- 
ment extends from 30 MeV to 100 GeV. It also gives a 
very good representation of the latitude distribution of 
the emission from the plane to the poles, and of the lon- 
gitude distribution. IC emission is a major component 
of the diffuse emission and dominates outside the Galac- 
tic plane and at energies below 100 MeV. The model 
reproduces simultaneously the 7-rays, synchrotron, CR 
secondary /primary ratios, antiprotons and positrons. In 
this sense it go es a long way toward s realizing our original 
goal, stated in IStrong et al.l l)2000(l . to reproduce astro- 
nomical and directly-measured data on cosmic rays in 
the context of a single model of the high-energy Galaxy. 

Obviously our optimized model is far from unique, 
both in the choice of parameters like halo size, diffu- 
sion coefficient, source gradient, cosmic-ray spectra etc. 
The purpose of this paper is to show that it is possible 
to construct at least one model which is consistent with 
all the relevant data within understandable limits on CR 
fluctuations and solar modulation. 

Based on the optimized model, a new EGRB spectrum 
has been derived (jStrong et al.ii2004y . 



We would like to particularly thank David Bertsch 
for assistance and discussions on the subject of the 
events and instrumental response of the EGRET tele- 
scope above 10 GeV and Seth Digel for providing the 
kinematically analysed H i and CO data used in this 
work. A part of this work has been done during a visit 
of Igor Moskalenko to the Max-Planck-Institut fiir ex- 
traterrestrische Physik in Garching; the warm hospital- 
ity and financial support of the Gamma Ray Group is 
gratefully acknowledged. Igor Moskalenko acknowledges 
partial support from a NASA Astrophysics Theory Pro- 
gram grant. Olaf Reimer acknowledges support from the 
BMBF through DLR grant QV0002. 



APPENDIX 
DESCRIPTION OF H i AND CO DATA 



The H I and CO data used in this work are based on more recent surveys than those used in IStrong fc Mattoa 
((1,996) and IStrong et alJ ((20001) . They were provided by S. Digel (private communication) who provided the following 
description. 

The annular maps are generat ed for 8 ranges of Galactocentric distance on the assumption of uniform circular 
rotation with the rotation curve of lClemena ((1985f) parameterized for Rq = 8.5 kpc, Vq = 220 km/s. Emission beyond 
the terminal velocity is assigned to the tangent point, and emission at slightly forbidden velocities in the outer Galaxy 
is assigned to the local annulus (7.5-9.5 kpc). The longitude ranges within 10° of ^ = and I = 180° are excluded 
from the integrations in all annuli. The boundaries of the of Galactocentric distance annuli are 1.5, 3.5, 5.5, 7.5, 9.5, 
11.5, 13.5, 15.5, 50 kpc. 

The CO data are from the 115 GHz line survey of iDame et alJ ((1987^ and the latitude range is |6| < 25°. The 
coverage is not complete within this latitude range but little or no significant CO emission is believed to be missed. 
The maps are of CO line intensity integrated over the (longitude-dependent) velocity range o f each an nulus and they 
have angular resolution 0.5° x 0.5° (set by the sampling pattern of the constituent surveys; see lDameet al. 1987). Th e 
units are velocity integrated radiation temperature (Wco), corrected to the intensity scale of lBronfman et all ((1988() . 
in K km/s. 

The H I data are a composite of several 21-cm line surveys (Table IX4]l. which were i nterpolated to a miiform grid. 
Calibrations were checked against the Bell Labs H i horn survey of St ark et all ((19921) . Brightness temperatures Tf, 
were converted to column densities of atomic hydrogen on the assumption Tgpin = 125 K uniformly. The few positions 
with Tb > 110 K had T^ truncated to 110 K. The maps have units of column density A^(II i)/10'^° atom cm"^. Angular 
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Table A4. HI surveys used for annular ring maps. 



Survey 




Angular resolution 
(HPBW), dog 


Region of Sky 


Weaver & Williams ("19731 
Heiles & Habine (1974T^ 


19791 


36' 
36' 

48' 

48' 




\b < 10°, I = 10-250° 
\b > 10°, <5 > -30° 




b < 10° , i = 240 - 350° 
b > 10°, 5 < -30° 



resolution is somewhat better than 1° and the maps extend to |6| = 40°. 

The Wco and iV(H i) maps described above were generated in 1996 using the then best available surveys of CO 
and H i. Since that tim e, surveys with gre ater coverage, angula r resolution, sensitivity, and improved calibration have 
been published, see Da me et al.1 l|2001j1 and lBurton et all l)1994j) . However these improvements would hardly affect the 
results presented in this paper. 
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